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Abstract 

We approach the study of non-integrable models of two-dimensional quantum field 
theory as perturbations of the integrable ones. By exploiting the knowledge of the 
exact S'-matrix and Form Factors of the integrable field theories we obtain the first 
order corrections to the mass ratios, the vacuum energy density and the S'-matrix of 
the non-integrable theories. As interesting applications of the formalism, we study 
the scaling region of the Ising model in an external magnetic field at T ~ T c and 



the scaling region around the minimal model M2J. For these models, a remarkable 
agreement is observed between the theoretical predictions and the data extracted 
by a numerical diagonalization of their Hamiltonian. 
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1 Introduction 



Two-dimensional integrable and relativistic models have provided a most valuable heuris- 
tic guide in the analysis of the infinite-dimensional space of Quantum Field Theories and 
the associated statistical models. The peculiar property of the dynamics that makes 
those models solvable consists in the existence of an infinite number of integrals of mo- 
tion. This feature manifests itself in several ways [|l], |2| . First of all, it leads to very severe 
selection rules: for instance, scattering processes accompanied by production or annihi- 
lation events cannot take place in such theories and therefore the only allowed scattering 
processes are purely elastic. Secondly, the n-particle scattering amplitudes can always be 
factorized into the product of the n(n — l)/2 two-particle elastic S-matrices. The two- 
particle S-matrix S ab (s) for the particles A a and A b , with mass m a and m b respectively, 
has only two elastic branch cut singularities at s = (m a — m b ) 2 and s = (m a + m b ) 2 
and fulfills by itself the unitarity equation S a b(s) S^(s) = 1 in the whole Riemannian 
surface of the Mandelstam variable s. Finally, the exact mass spectrum and the com- 
plete set of S-matrices of integrable QFT may be computed according to the so-called 
bootstrap principle which states that the bound states have to be regarded on the same 
footing as the asymptotic states. The drastic simplification of the on-shell properties 
of the integrable field theories has far-reaching consequences since it greatly facilitates 
both the computation of their thermodynamic quantities and the determination of their 
correlation functions (see, for instance [10-20]). For instance, the usual lengthy compu- 
tation of matrix elements of local operators which enters the spectral representation of 
correlation functions can be shortened in the case of integrable models to solve a finite 
number of functional and recursive equations. At present, the exact S'-matrix, the exact 
mass spectrum and the correlation functions of a large number of integrable models have 
been determined^ and confirmed either by analytic lattice computations or by numerical 
simulations [32-39]. One of the most remarkable examples of statistical model solved by 
using the bootstrap approach is represented by the Ising model in a magnetic field at 
T = T c : the exact S'-matrix proposed by Zamolodchikov has been the starting point 
-^or a brief list of references see [1-20]. 
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for the solution of the long-standing problem of determining the spin-spin correlation 



function of this model [11 



Despite all the distinguishing qualities of integrable models and their successful ap- 
plications to statistical mechanics, the actual realm of particle physics is however that 
of non-integrable QFT. These theories generally present the striking phenomenon of res- 
onances as well as multi-production reactions. Many interesting statistical models in 
the vicinity of their fixed points may be described by non-integrable QFT and, to make 
any progress in the computation of their thermodynamical quantities (or, simply to un- 
derstand better their qualitative features), it would be highly desirable to develop an 
appropriate formalism to deal with the lack of integrability. This task is notoriously dif- 
ficult because the rich physical scenario of the multichannel systems, represented by the 
non-integrable field theories, is usually accompanied by great mathematical complexities. 
In fact, once one has given up the integrability condition, the infinite number of thresh- 
olds associated to the production processes greatly influences the analytic structure of 
the scattering amplitudes, inducing a rich pattern of branch cut singularities in addition 
to the pole structure generally associated to the bound or resonance states (Figure 1). 

The mathematical difficulties of the problem are well known and cannot be easily 
circumvented, as illustrated by the following line of reasoning. Let T\ < T 2 . . . < T n . . . 
be the infinite sequence of energy thresholds of the scattering amplitudes entering the 
S'-matrix. In each interval T n < s < T n+1 , the physical S-matrix is given by an n x n 
unitary matrix S^l (s) = 1, . . . , n), where here the suffixes % and j are multi-particle 
collective indices that generally label the different channels of the scattering process. 
At each threshold, a new channel opens up and correspondingly, more matrix elements 
have to be introduced. Due to the nested structure of the S'-matrix as a function of 
s, i.e. S (1) -> S {2) ... -> S (n) , one may be inclined to approach the difficult problem of 
computing the infinite number of scattering amplitudes by applying a recursive procedure. 
Namely, assuming that all the n 2 matrix elements S^\s) have already been determined 
as analytic functions of s in the interval T n < s < T n+1 , one might try to analytically 
continue them into the next interval T n+ \ < s < T n+ 2 and then, as consequence, compute 
the new unknown scattering amplitudes entering the larger matrix by using the 
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unitarity equations relative to the new range of energies. There is, of course, an important 
flaw in the above argument: the recursive procedure is in fact based on the assumption 
that the matrix elements S^\s) can be prolongated through the different thresholds 
without being halted by the presence of singularities, an assumption which is obviously 
false. Hence, the complexity of non-integrable QFT consists in the fact that in any 
energy interval one must take simultaneously into account all the inelastic thresholds of 
the theory. This feature is explicitly shown by a simple example discussed in Appendix 
A. For a discussion on the inter-related aspects of analyticity and integrability, see also 



The problem of analysing non-integrable field theories (or at least a class thereof) is 
not however as impracticable as it may seem at first sight. A possible approach to their 
study is suggested by the observation that although integrable and non-integrable QFT 
possess quite different properties as far as scattering amplitudes and mass spectrum 
are concerned, their ultraviolet limit may be described by the same Conformal Field 
Theory (CFT). From this point of view, an integrable field theory is nothing more than 
a particular deformation of a Conformal Field Theory whose euclidean action may be 
written as 



where is one of the relevant operators which mantains the original integrability of 



two ways of defining a non-integrable field theory: (a) the first possibility is to simply 

add a perturbation induced by a non-integrable operator to the CFT action ; (b) an 

alternative possibility is to deform the CFT action by means of a linear combination 

of several operators in which at least one of these or all of them are individually 
2 Zamolodchikov has outlined some criteria on the integrability of the deformations of CFT. For 
a generic minimal model, for instance, integrable QFT are obtained as deformations by the relevant 
operators $1,3, $1,2 and $2,1 [||. 



ref . |7| . 





the conformal mo delft Starting now from a given conformal model, there are at least 
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integrableQ. For the corresponding action in the last case we have 

A = Acft + T,9j J d 2 x<S> 3 (x) . (1.2) 



Apart from numerical investigations or perhaps standard analyses based on con- 



formal perturbation methods, non-integrable QFT of type (a) are presently difficult to 
study and they will not be considered in this paper. On the contrary, non-integrable 
QFT of type (b) will be shown in this paper to be suitable of a purpose of theoretical 
analysis. This is particularly important because they provide in many respects the most 
interesting realization of non-integrable models. Let us illustrate the key observation that 
permits their analytic approach. 

In order to study the non-integrable QFT associated to the euclidean action (|1.2|) , it 
is convenient to regard the latter as the action relative to a deformation of an integrable 
QFT rather than a multiple deformation of CFT. This consists in grouping differently 
the terms in ( |1.2|) and rewriting it as 

A = A* + Y,9iJ d 2 x^(x) . (1.3) 

There are several advantages to adopting this point of view. 

• Firstly, going to the Minkowski space, we can start our analysis of non-integrable 
QFT by using the particle basis of the integrable model associated to A\ nt . Although 
their mass spectrum will not generally coincide, the particle basis of the latter model 
is surely more suitable than the conformal basis for capturing the large distance 
properties of the non-integrable theory^. 

• Secondly, the integrable QFT with which we start our analysis are solvable theories, 
although not necessarly free. This means that we know how to compute the exact 



3 In the case when all the operators $i are individually integrable but their simultaneous presence 
renders the field theory non-integrable, it can be argued that this is due to different null-vector structure 
of the fields $i Q or it can also be inferred by numerical studies of the spectrum, as discussed in the 
next sections. 

4 This is generally the case for massive field theories. For massless non-integrable field theories the 
analysis may be more subtle and will not be pursued here. 
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matrix elements of all their local operators, in particular those entering the action 
Q1.3Q . Hence, we can develop in this case a perturbative approach to non-integrable 
QFT based on the exact Form Factors of the integrable ones. As we will discuss 
in this paper, this approach seems generally more powerful and efficient than the 
one based on usual conformal perturbation theory. In fact, a large set of quan- 
titative information -concerning the mass corrections, the renormalization of the 
vacuum energy density and the S'-matrix of the non-integrable models- can already 
be obtained at the lowest order in the coupling constant gj with a high level of 
reliability^. 

• Finally, in the case of multiple deformation of CFT by separately integrable oper- 
ators, we can choose any of them for defining our initial solvable QFT. Obviously, 
each choice leads to a selection of particle basis, bound state structure and dy- 
namics, on which we can construct our perturbative approach. However, since the 
actual dynamics of non-integrable QFT must be independent of the solvable model 
we start with, this implies that perturbative expansions based on different Form 
Factors and spectra must be related each other by a set of mathematical identities. 

The most natural interpretation of eq. (|1.2|) when all are individually integrable 

operators, is then that of an interpolating action between exactly solvable models, the 
interpolation being realized by varying the relative values of the coupling constants 
(Figure 2). A remarkable example of such interpolating theory is provided by the Ising 
model away from the critical point, i.e. at T ^ T c and in a presence of a magnetic 
field h. In the pure magnetic direction, the spectrum of this model consists of eight 
stable particles, five of which above threshold ||. In the pure temperature direction, 
on the contrary, the spectrum presents only one massive state, which can be interpreted 
either as a kink (in the low-temperature phase) or as a particle (in the high-temperature 
phase) [22-26]. The interpolation between the Hilbert spaces constructed on these rather 
different bases and the rich physical scenario associated to it will be analysed in Section 
5. 

5 The formulas which we will derive in Section 3 can be regarded as "Born approximation formulas" 
for a non-integrable QFT. 
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At this stage of the discussion, the careful reader may have noticed that we could have 
directly introduced the non-integrable QFT by means of the action ( |1.3|) without further 
specifying how the integrable action A\ at has been actually obtained. This is indeed true 
and in fact, all the formalism set up in this paper to study non-integrable models basically 
relies on the properties of relativistic integrable models as defined, for instance, by the 
framework of the analytic theory of the S'-matrix. Although the perturbed CFT approach 
we have used in this section may then appear only as a convenient way to introduce and 
define non-integrable QFT, however the eventual knowledge of the underlying conformal 
model actually adds an important piece of information on the operator content and the 
ultraviolet properties of the theory. In fact, by the knowledge of the CFT ruling the 
short-distance Operator Product Expansion of the theory, one can severerly constrain 
its high-energy behaviour and obtain, among other things, a bound on the cross-sections 
of non-integrable QFT. This aspect as well as others of non-integrable models will be 
discussed in a separate publication . 

The paper is organized as follows. In Section 2 we define the perturbative scheme 
based on the so-called intermediate state representation and the Form Factors of inte- 
grable theories. We also discuss the necessity of introducing some counterterms in the 
perturbative series in order to have a consistent formulation of the theory. In Section 3 we 
specify our analysis to the case of two-dimensional integrable QFT and derive the main 
formulas which rule the variation of the mass ratios, the vacuum energy density and the 
S'-matrix. Section 4 is devoted to the simplest model which can be discussed by means of 
these techniques, namely the minimal model A^2,7 perturbed by the operators $ 1>2 and 
$1,3. Another example of non-integrable quantum field theory is given by the Ising model 
away from criticality, at a generic point in the bidimensional phase diagram spanned by 
the magnetic field h and the reduced temperature (T — T c )/T c . In view of its important 
role in statistical mechanics, the Ising model away from criticality will be discussed in 
some detail in Section 5. In particular, we will compute the first order corrections to the 
spectrum of the theory when a small thermal perturbation is added to the pure mag- 
netic case. Our conclusions are then reported in Section 6. Three appendices complete 
the paper: the first one shows how the inelastic thresholds enter the elastic scattering 
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amplitude of a simple example of quantum mechanics; the second appendix discusses 
the Faddeev-Zamolodchikov algebra and the disconnected terms present in the matrix 
elements of local operators; the last one gathers the relevant Form Factor expressions for 
the $1^ deformation of the M( 2 j) model. 

2 The Intermediate State Representation 

Let us consider a quantum field theory in the Minkowski space, defined by the action 

A = Ao + Ai = A -X J d 2 x *(ar) , (2.4) 

where Ao denotes here the minkowskian action of the unperturbed theory and ^ one of 
its operators. We suppose that the QFT associated to the action Ao is exactly solvable 
(although not necessarily free), i.e. we assume that the spectrum of particles, their 
scattering amplitudes and the matrix elements of the operators of the theory (and in 
particular those of ty) are all known. For the sake of simplicity, we consider in this 
section the case of isospectral perturbations of a solvable theory. This means that the 
spectrum of the theory described by the total action A will be made of the same number 
of particles of the unperturbed one; in other words, the new interaction Ai only effects 
the values of the masses of the physical particles but not their stability properties^. 

Let us now describe the properties of the theory associated to the action A. Under 
the hypothesys that the interaction term is turned off at t — > ±oo, it is possible to adopt 
the formalism of the asymptotic "in" and u ouf states. We are interested in computing 
the scattering amplitude 

S{qi, ■ ■ ■ , q n - q'i, ■ ■ • , q'J = • • • > <4l?i, • • • > Qn) m = in (q[, • • • , ^ m \S\q h . . . , q n ) in , 

(2.5) 

where qi and qtj label the momenta of the in-going and out-going set of particles. Since 
in the remote past t — > — oo the interaction is not present yet, the asymptotic in states 
coincide with the unperturbed ones 

, qn) m = \qi, ■ ■ ■ , qn) n ■ (2.6) 

6 We will comment on the more general case in Section 5 when we will consider the Ising model. 
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As usual, the scattering operator S appearing in (|2.5| ) can be obtained as the limit 

S= lim U(t,-t) (2.7) 

of the time evolution operator U(t,t ), which is the solution of the equation 

ijU(t, t ) = H U{t, t ) , U(t , to) = 1 , (2.8) 

where H = Ho + Hi denotes the Hamiltonian associated to the theory fl2.4|) . Following 
the standard quantum mechanical procedure (see, for instance |[21|), the operator U can 
be factorised as U — U Ui, where U and Ui are solutions of eq. (|2.8j ) with H replaced 
by Hq and Hi(t) = Uq 1 HiUo ) respectively. Then, we can write the scattering operator 
of the theory ( p^) as S = SoSi, where So = lim^+oo Uo(t, —t) is the unperturbed and 
exactly known scattering matrix 

So{qi, ■ ■ ■ , Qn Ql, ■ ■ ■ , q' m } — ^oWn ■ ■ ■ > ^ml^b • • • ) In)^ = ^(Ol) • • • j Qm\^o\<ll, ■ ■ ■ , 9n)o" J 

(2.9) 

while Si has the usual formal representation 

Si= lim Ui(t,-t) =Texp(iA z [*]) . (2.10) 
The scattering amplitude is therefore given by 

° ut (q[, ...,q'Jqi,..., q n ) m = out (q[, gjTexp (-zA | rf 2 xvl/(x)) . . . , = 
= E L TT L / ■ ■ ■ • • • • sJ»l T • • • I*. • • • ^)o" . (2- 11 ) 

fc=0 

where ( |2.9| ) has been used in order to absorb the factor Sq. 

In ordinary lagrangian perturbation theory based on free theories, the computation of 
scattering amplitudes would now proceed through the use of creation and annihilation op- 
erators and Wick's theorem, finally leading to the diagrammatic expansion characteristic 
of Feynman covariant perturbation theory. Although this approach cannot be generally 
followed here (because we might not know, in general, a local lagrangian formulation 
of the theory associated to ^4o) ; the knowledge of the exact solution of its dynamics, 
in the form specified at the beginning of this section, naturally suggests to compute 
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the scattering amplitudes ( |2.11| ) within a framework quite analogous to that of ordinary 
time-dependent perturbation theory in quantum mechanics. In other words, we initially 
insert between the operators ^(xi) and ^(xi + ±) (I = 1, . . . , k — 1) in the second line of 
( |2.11| ) a sum over a complete set of asymptotic states of the unperturbed theory 

£ |n)S* * <n| = 1 = 5»8* «*<n| , (2.12) 

n n 

with \n) denoting an asymptotic state containing n on-shell particles. Then, the inte- 
grations on the space coordinates in ( |2.11| ) are immediately performed and lead to delta 



functions constraining the total momentum of the intermediate states to coincide with 
that of the initial and final states. In doing the integrations on the time variables, the 
time ordering prescription gives rise in this case to the appearence of energy denominators 
so that one ends up with the following expression 

° ut (q[, ...,q'Jqx,---, q n ) m = ^M, . . . , <4|<2i, . . . , q n )t + 

(m n \ 

E^-E?i {-a ou M,---,q , m \m\qi,---,qn)t l + (2.13) 
3=1 i=i / 



X 



2™f2 i^"'r^A( E - E l+ ie )---( E - E k-l+ ie 

x ^ ll ...,g^|*(0)|Th>o...o<"»-i|*(0)|?i,...,?„>J 1 ]} , 

where E and Ei (Q and Pi) denote the total energy (momentum) of the initial state and 
of the 2-th intermediate state, respectively, and each intermediate sum can be equiva- 
lently taken either on the basis of the in states or on that of out states. Since the matrix 
elements between asymptotic states of the perturbing operator ty(x) are supposed to be 
known, the scattering amplitudes in the quantum field theory ( |2.4j ) are in principle com- 



putable quantities, order by order in the coupling constant A. The above expansion over 
intermediate states must be contrasted with the usual formalism of covariant perturba- 
tion theory in which both energy and momentum are conserved in the internal lines of 
Feynman diagrams where, however, such lines correspond to off-shell particles. 

As it is, the above formula (|2.13 ) is not completely correct though. In fact, the new 



interaction changes both the vacuum energy density and the mass of the particles. Hence, 
we have to refine the action Ai by introducing some counterterms to take properly into 
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account the correct normalisation of the states. We require the validity of the following 
normalisation conditions for any value of the coupling constant: the normalisation of the 
vacuum state 

(0|0> = o<0|0) = 1 , (2.14) 
and the normalisation of the one-particle states 

out {g , lq) in = out^^n = ^ gtf _ g l) . (2 .15) 

The two above conditions should be enforced order by order in perturbation theory when 
using ( |2.13| ) to compute the vacuum to vacuum and the one-particle amplitudes. 

Equation (|2.14j) leads to subtract a constant term 5£ vac (X) from the interaction den- 



sity. This extra term obviously measures the variation of the vacuum energy density 
under the effect of the perturbation. This variation is usually ignored in lagrangian per- 
turbation theory where the prescription of disregarding the disconnected vacuum bubble 
diagrams is adopted. We keep track of it here because it is a measurable quantity for the 
class of models we are considering in this paper. 

The correct one-particle normalisation may be obtained by introducing a 'mass' term 
operator in the interaction density. This operator, denoted here by O^(x), can be defined 
in terms of its (unperturbed) Form Factors, given by 

F^ = {0\O^(0)\q 1 ,...,q n )t l = S n>2 . (2.16) 

With this definition, the coefficient in front of the operator 0^ 2 \x) in the interaction 
density has to be determined by imposing eq. Q2.15 ) order by order in the coupling A and 
plays the role of a mass counterterm 5m 2 (A). 

In summary, the correct formula for the scattering amplitude is given by 

° ut Wn ■ ■ ■ , Q'JQU ■ ■ ■ , q n ) m = ^M, ....<L<1. q n )P + 

(m n 
3=1 3=1 

X { . . . , q'J (A* (0) - l -5m 2 0^(Q) - SsJ) \q x , . • • , q n )™+ (2.17) 

+ E le-^+i) ^ ■■- q,J " \ 5m2 ° m ^ ~ ^vac) |ni> X 

(n 1 \^(0)-^5m 2 O^\0)-S£ vac ^\q 1 ,...,q n ) l n + ...} . 
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Let us remark that the above expansion appears as the most physical one since it deals 
with the true physical degrees of freedom of the problem. However, as usual in quantum 
field theory appearences of divergent contributions are expected when the above formula 
is applied beyond the first perturbative order. A general discussion of such divergences 
and of the renormalisation procedure which must be adopted to deal with the infinities, 
seems to be an interesting open problem in the unconventional setting we are considering. 
In fact, it must be taken into account that for the perturbations of the interacting theories 
we are interested in, the perturbing operator has in general non- vanishing matrix 
elements on all the asymptotic states. This means that a resummation of an infinite 
number of terms is required at any perturbative order beyond the first. The general 
aspects of this problem will not however be further investigated in this paper. Rather, 
we will concentrate our attention here on the first order approximation in order to study 
the effects induced on an integrable model of two-dimensional quantum field theory by 
a small perturbation which breaks the integrability. 

3 Perturbing Two-dimensional Integrable Models 

In this section we apply the perturbative scheme exposed in the previous section to the 
case in which the term Aq in the euclidean version of the action (|2.4j) corresponds to 
an integrable deformation of a two-dimensional conformal field theory through a relevant 
operator $(a;), i.e. 



Before proceeding in our analysis, a brief review of the main properties of the integrable 
relativistic field theories on which we base our sequent considerations seems in order. 

3.1 Basic Features of Integrable Theories 

Integrable theories are characterized by the existence of an infinite number of integrals 
of motion. This circumstance gives the possibility of solving them along completely non- 
perturbative methods. The first consequence of the infinite number of conserved currents 
is that the scattering processes in the integrable theory are completely elastic, i.e. the 




(3.1) 
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final state contains the same number of particles with the same momenta of the initial 
one: the only results of the scattering processes are then a possible exchange of the 
"color" quantum numbers between particles in the same mass multiplet together with 
the elastic phase shifts associated (semi-classically) to the time delays of the wave packets 
[|l| . Moreover, in these cases the S-matrix is completely factorised into the product of the 
two particle scattering amplitudes defined by 

mit (c(9 3 )d(9 4 )\a(9 1 )b{9 2 ))^ = (2vr) 2 5(^ - 9 3 )6(9 2 - 04)^1(01 - e 2 ) , (3.2) 

where the rapidity variable 9{ parameterises the on-mass shell energy and momentum of 
the particles as (p°,p[) = (m ; cosh 9 iy rrii sinh #j). Since the Mandelstam variable s, given 
by 

s ab {0i - 6 2 ) = {p a (9 x ) + p b {6 2 )) 2 = m 2 a + m 2 b + 2m a m b cosh(^ - 9 2 ) , (3.3) 

is the only independent relativistic invariant in the process, the scattering amplitude 
depends on the rapidity difference only. Once combined with the general principles of 
analyticity, unitarity and crossing symmetry, the elasticity and factorisation properties 
usually lead to the exact determination of the S-matrix in integrable models (see for 
instance [1-9]). 

The second important characteristic of the two-dimensional relativistic integrable 
models is that, in addition to the S-matrix, we can also determine exactly the matrix 
elements of the local operators 0(x) of the theory on the asymptotic states, i.e. 

^.. bm FZ.. an (e[, . . . , e'je,, . . . , e n ) = . . . uOWM^i) • • • <o n ))* . 

(3.4) 

We can always restrict our attention to those matrix elements with no particle on the 
left side, the so-called Form Factors (FF) 

FZ..a n (0n . • • A) = 0(010(0)1^(0!) . . . , a n (6 n )y o n . (3.5) 

In fact, the generic matrix element ( |3.4j ) can be obtained in terms of the Form Factors 
through the analytic continuation 

bi...bm F ai...a n (^l: ■ ■ ■ > Q'mftu • • • > On) = F b\...b mai ...a n (0[ + ilT, . . . , 9' m + in, 9 1: . . . , 9 n ) 

+ disconnected parts , (3.6) 
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where the bar denotes the charge-conjugated particles. The "disconnected parts" in the 
right hand side of (|3.6| ) obviously appear when some of the primed rapidities coincide 
with the unprimed ones (this point is discussed in more detail in appendix B). 

As well known, the computation of the Form Factors can be performed once the exact 
S'-matrix and the bound state structure of the theory are known. In fact, the monodromy 
properties of the FF are constrained by the so-called Watson equations, given by 

F ai...a i a i+1 ...a n {^U ■ ■ ■ > #i> ^i+l> • • • > ^n) = 

= S am+\(0i - F Z-b l +ib l ...a n ( e ^ " " " ' ^+1> &h ■ ■ ■ , #n) , 

Ka^x + 27rz, e 2 ,..., e n ) = FZ... anai (e 2 , ...,e n , e x ) , (3.7) 

whereas their analytic structure is closely related to the underlying pattern of singularities 
induced by the (multi)-scattering processes |TT|, [12], |13| and is ruled by a set of recursive 



equations. The simplest of these equations^ are obtained by the residue of the FF at the 



kinematical poles of particle-antiparticle singularities of relative rapidity 6 = in [|13 



- i lim Fg^ (6 + 6,6,,..., 6 n ) = ( 1 - J} S aai (6 - 0<) J F^.. an (6 U . . . ,6 n ) . (3.8) 

Additional recursive equations are induced both by the bound-state poles and higher- 
order poles due to multi-scattering virtual processes [11], [12], [U| . Finally, the asymptotic 
behaviour of the FF for very large values of the rapidities can be controlled according to 
a simple criterion proposed in 0] : denoted by xq the scaling dimension of the operator 
0(x) and by yo the real quantity defined by 



a I ~, 11— Or, 



we have 



vo<y ■ ^ 

The equations and the constraints which we have briefly illustrated above prove in 
general sufficient to determine the FF of the local operators of the integrable theories 
(see, for instance [11-20]). 

7 For simplicity we quote the formula corresponding to the case of diagonal scattering. 
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3.2 First Order Perturbation Theory 

Let us now deform the integrable action (|3.1|) by adding to it the relevant operator $?(x). 
Both the operators <$>(x) and ty(x) are assumed to be scalar and relevant, with their 
scaling dimensions denoted by x$ and x^, respectively The theory (|2.4j) depends in 
this case on the two dimensionful couplings constants^ g and A. Since g ~ M 2 ~ Xi ' and 
A ~ M 2-2 "*, where M is a mass scale, we can decide to use g as dimensionful parameter 
of the theory and the dimensionless combination 

X = \g- 2 ^ (3.10) 

as a label of the different Renormalization Group trajectories which originate from the 
fixed point at g = A = (see Fig. 2). For example, if N(x) denotes the number of stable 
particles in the spectrum of the theory, their masses may be expressed as 

m a (g, X ) = C a (x)g^ , a = l,2,...,N( X ) , (3.11) 

where C a (x) is an amplitude which characterises the whole trajectory. Similarly, the 
vacuum energy density can be written as 

S vac (g,x)=V(x)9 T * 5 ■ (3-12) 

Dimensionless quantities, as for instance mass ratios, only depend on x an d therefore 
they do not vary along the trajectories of the Renormalization Group. 

Once the new interaction A / d 2 x^(x) is switched on in the action, the integrabil- 
ity characterising the unperturbed theory is generally lost and the S-matrix becomes 
extremely more complicated. Inelastic processes of particle production are no longer for- 
bidden and the analytic structure of the scattering amplitudes will present additional cuts 
due to the higher thresholds. In particular, their expression is no longer factorized into 
the sequence of two-body scattering amplitudes and, even in elastic channels, the only 
surviving restriction on the final momenta comes from energy-momentum conservation. 

8 For the sake of simplicity of notation, we assume that no other interaction term is generated by 
renormalisation effects (this is the case of the two models which we will explicitely discuss in this paper) , 
although the first order results of this section do not actually depend on this assumption. 
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The knowledge of the matrix elements ( |3.4| ) of the perturbing field ^>(x) ensures the 
possibility to compute perturbatively the amplitudes of the inelastic processes as well as 
the corrections to the elastic ones. To the first order in A and with an obvious extension 



of the notation, equation (|2.17|) reads 



°*(bi(d) ■ ■ ■ b m (q l m )\aM) ■ ■ ■ «n(^)) m ^ (3.13) 

* Smn ° Ut o(bM) • • • b^WM) ■ ■ ■ <PI))T + 

(n m \ 

EK-E9fc "ofafoJ) ■ ■ ■ (*(«)+ (3-14) 

k=l fc=l / 

i N \ 

- - E 6M 2 ab O%j(0) - ^ vac | 0l (rf) . ..aM))o ■ 

Z a,6=l / 

The "mass operator" O ab (x) is defined assigning its form factors, which with an obvious 
generalization of eq. ( P-16[ ), are given by 

(2) 

S ba2 ■ (3.15) 

The first order corrections to the masses of the particles and to the vacuum energy density 
are then obtained imposing the conditions ( |2.15| ) and (|2.14 ). As a result, we have 



5M-1 ~ 2\F- b l(i>K,Q)5 mamb , (3.16) 

SE^-X [o(0|*|0) ] . (3.17) 
Care must be however taken when using the rapidity parameterisation in eq. (p.!4|) . 



Let's illustrate this point by considering the first order correction to some elastic process 
ab — > cd. In the unperturbed theory, this process is characterised by the scattering 
amplitude S^f(8), where 9 = Q\ — 6 2 denotes the rapidity difference of the colliding 
particles. In two dimensions, the momenta of the particles in a two-body elastic collision 
are individually conserved even in absence of integrability, so that the general elastic 
amplitude S^ b (9, x) can be introduced through the relation 



out (c(9 1 )d(e 2 )\a(9 3 )b(9 4 )) m = (2tt) 2 5(9 1 - 9 3 )5(9 2 - 4 ) Sf h {9 x - 9 2 , x ) ■ (3.18) 

Notice, however, that away from the integrable direction (i.e. \ = 0), the scattering 
amplitude S^b{9, x) is no longer a meromorphic function of 9 since the opening of inelastic 
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channels induces an additional analytic structure. When we compute the correction 
to S^f(9, x) around x = 0, we must take into account that, since the total energy of 
the colliding system is fixed, the variation in the masses given by eq. ( 3.16|) induces a 
corresponding change in the rapidity difference expressed by 

m a 5m a + m b 5m b + {m b 5m a + m a 5m b ) cosh 9 
69 " m a m b sinh6 " (3 ' 19) 



Then the correction to the amplitude can be decomposed as 

S X ■ (3.20) 



89 d\ 

x=o 

The first order result for this quantity is obtained by using formula Q3.14Q . Taking into 
account the cancellation occurring between the disconnected parts in eq. (|3.6|) and the 
contributions of the counterterms, one finally obtains 

5S? b (9, X )~-i\ F ^ ab{6 \ ft , (3.21) 
m a m b sinn 9 

where 

*£*(0i - 2 ) = F^ ab {9, + in, 9 2 + in, U 9 2 ) . (3.22) 

The right hand side of ( |3.21| ) employes the expression of the Form Factor at very special 
values of the rapidity variables. According to eq. ( |3.8| ), the Form Factors present pole 
singularities whenever the rapidities of a particle-antiparticle pair differ by m and, in 
fact, these kinematical poles are often explicitly inserted into the denominator of their 
parameterization. Apart from a term encoding the monodromy properties, this parame- 
terization may be written as Q/T> where both Q and D are polynomials in the variables 
cosh 6^-, the denominator being fixed by the pole structure whereas the numerator de- 
termined by means of the residue equations, as for instance those of eq. ( p.8|) . From the 



finiteness of the left hand side of eq. ( |3.21| ), we expect therefore that the "m singulari- 
ties" of the denominator of the Form Factors F3 J9±, 9 2 , 9 3 , 9$) should be cancelled by 
the polynomial Q, once evaluated at the specific rapidity configuration of eq. Q3.22p . This 
prediction will be explicitly checked in the next section and should hold in general when- 
ever the perturbing operator is local with respect to the fields which create the particles 
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in the unperturbed theory^. 

Eqs. (|3.16|) ) (|3.17|) and Q3.21|) are the main results of this section. The best use of 
these formulas is to get rid of the explicit dependence on the normalisation of the per- 
turbing operator by defining universal quantities, as for instance ratios of the mass shifts. 
Hence, under the validity of the linear approximation, all the universal quantities of non- 
integrable field theories can be entirely expressed in terms of the Form Factors of the 
integrable ones. Comparison of the theoretical predictions with their numerical determi- 
nations will be presented in the next sections of the paper. 

It is particularly instructive to specialise the above discussion to the "trivial" case 
in which the perturbing operator coincides with the operator that defines 

the initial integrable theory. In this case of course the physics should be invariant since 
the result of the additional perturbation simply corresponds to a shift of the coupling 
constant of the original integrable model by an amount 5g = A. The variations of the 
masses of the particles and of the vacuum energy density corresponding to such a shift can 
be directly computed from eqs. (|3.11|) and ( |3.12| ), respectively. On the other hand, we can 
also apply our general formulas ( |3.16| ) and ( |3.17D to estimate the first order corrections. 
The two different routes coincide as long as the following identities are valid 

aaK J a (3.23) 

£vac = ^(0|e|0> , 

where Q(x) = 2ixg(2 — x$) is the trace of the energy- momentum tensor for the 

trajectory x = 0. The above two relationships are indeed true and can be easily derived 
by other means. However, it is interesting to notice that their validity emerges in this 
context for the role of consistency equations which they play. Obviously, by the same 
token we can generate an infinite number of identities involving Form Factors of the 
original integrable field theory by considering higher multi-particle scattering processes. 
For instance, next to Q3.23 ), a new identity is obtained by comparing eq. ( |3.2U ) with 
eq. fl3.21| ): since x is constant in the case we are considering, we have 

dsim _ 1 f%jo) 



86 2ni s ab {6) 



(3.24) 



9 We defer to the section devoted to the Ising model the discussion of the case in which this condition 
is not fulfilled. 
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This identity provides a simple and unique way to normalise the four-particle Form 
Factors of the stress-energy tensor. It may be then particularly useful in the study of 
massless field theories where the first relationship in eq. ( |3.23| ) cannot be used for this 
purpose. For instance, it is easy to check that eq. fl3.24|) applies to the Form Factor of 
the massless model considered in ref. 

It is also obvious that the first order inelastic amplitudes computable by formula 
( |3.14| ) must vanish identically when we choose ^>(x) = This is ensured by the 

fact that the form factors of the stress-energy tensor F® a „(0i, • • • , n ) factorise the term 
P^P^, with P^ = Yn=iPi denoting the total energy-momentum of the set of particles. 
Since pf — *> — p% when the i— th particle is crossed from the initial to the final state, P M is 
zero for a set of particles entering a physical scattering process. Only in the case of elastic 
scattering, the zeros coming from the factor P^P^ are cancelled by the kinematical poles 
and therefore relations analogous to eq. ( |3.24j ) are obtained. 



4 Non-Integrable Deformations of the Minimal Model 

M m 

As already indicated in the introduction, the main idea of using a perturbative expansion 
based on the Form Factors of integrable theories is because we expect that this kind of 
series should be capable of approximating the dynamics of the non-integrable field the- 
ories close to the integrable ones better than any other approach. One of the physical 
reasons for this expectation is that the integrable field theories should provide from the 
start the right multi-particle basis in the Hilbert space of the perturbed, non-integrable 
ones: hence, the differences between their physical properties are presumed to be small 
and calculable. Said in mathematical terms, such perturbative expansion should be par- 
ticularly significant since the exact expression of the Form Factor already corresponds to 
a resummation of an infinite number of terms originating from the action ( |1.2| ). However, 
the final conclusion on the efficiency of the Form Factor perturbative approach should 
come from some explicit and direct comparison with a set of data obtained from other 
sources. This is what we are going to accomplish in this section as well as in the next one. 
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As our first example, we will discuss the scaling region around the fixed point described 
by the non-unitarity minimal model Mt^j)- This model is particularly appealing for its 
simplified dynamics whereby the significant physical effects we are looking for will not be 
masked by other additional complications. 

The minimal conformal model M( 2 j) has only two primary operators, $13 and $1,2, 
both of them relevant with scaling dimensions given by —6/7 and —4/7 respectively 
[p8| . The perturbations of the conformal action either by the operator $i ; 3 or by the 
operator $ 12 are both known to be, separately, integrable [0, f|. In their massive phase, 
both perturbations have two stable massive particles denoted by Ai{9) and A 2 (6), with 
a mass ratio and a scattering matrix which depend, however, on the integrable direction 
considered. The exact two-particle elastic S-matrix and other relevant information about 
the two different integrable deformations are summarized in the following table 



$1 3 deformation 



Sum 

S 12 (9) 
S 22 (») 



(I) 



(I) 

m 



mi = 2 cosf = 1-6180. 

mi 5 



P 



-0.1314..m? 



$1 2 deformation 



3 n (0) 



(1) 



1112 

mi 



2 cos^ = i- 9696 - 



8(sin|+sinf-sin^j 



-0.2211. .ml 



(4.1) 



where 



la) = 



tanh \ {Q + ma) 
tanh \ {9 — ma) 



and the bound state poles in the S'-matrix amplitudes related to the particles A{ are 
identified by the index i placed above the functions (a). 

At the present, the $1 3 perturbation is the most studied of the two and in particular 



the Form Factors of the two primary fields are known flTJJ. This suggests to adopt as our 
initial euclidean action Aq that relative to the field $1,3, i.e. 



A> = -4(2,7) + g J d 2 x $1,3(3;) 



(4.2) 
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and then add to it the other deformation Q 

A = A + A J d 2 x $ lt2 (x) . (4.3) 



It is natural to assume that the spectrum associated to the action (|4.3|) consists for all 
positive values of the dimensionless parameter x of two massive non-degenerate excita- 
tions (the validity of this assumption can be directly confirmed by means of a numerical 
technique that we will mention below). Hence, the QFT described by ( |4.3|) is of the type 
of isospectral theory discussed in Section 2. This circumstance, together with the rela- 
tive simplicity of the scattering theory in both integrable directions, makes this theory 
the ideal playground for testing the Form Factor perturbative scheme. In order to make 
explicit predictions, we have listed in Appendix C the matrix elements of the operator 
$12 that we need in the sequel. 

Let us initially compute the first order correction to the elastic scattering amplitude 
Sn(6). According to the general formula ( |3.21| ), this is given in terms of the four- 
particle Form Factor given by eq. (|C.4|) , computed at specific values of the rapidities. The 
corresponding particle configuration consists of neutral particles with rapidities differing 
by in, i.e. the FF is on the resonant configuration ruled by the "kinematical poles" of 
the factors cosh(9ki/2) in the denominator of ( |C.4j) . However, the factor Q(%) present 
in the numerator of the four-particle Form Factor exactly cancels these divergences and 
produces finite result 

tanh±(# + ^) Ac ( nn . 2 tt\ cosh 2 § 1 + 2 cos f cosh 9 



5ll( ^ A) ^tanhi(^-zfj VM) 2 l 32sin sj sinhtf (sinh 9 - i sin f )' " (4 " 4) 

The rapidities employed in the above formula are obviously defined on the unperturbed 
mass-shell condition. The finiteness of the above formula is an explicit demonstration 
of the cancellation of the kinematical singularities of the FF discussed in the previous 
section. 



10 Throughout this section we will refer to the region of the coupling space of the theory (4.3) in which, 
within the standard CFT normalisation of the fields, g is real and negative and A is purely imaginary 
with positive imaginary part. The presence of the imaginary coupling is not surprising in view of the 
non-unitary nature of the original conformal model. 
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The correction term in ( |4.4| ) exhibits a second order pole at 9 = 27r/5. The increased 
order of the pole signals that the position of the bound state pole in the #-plane has 
become in the non-integrable theory a function of the coupling constant A. This suggests 
the reabsorbing of the first order term in ( |4.4| ) as a correction to the location of the 
original simple pole, given now by 

2n Xcq 27T 7T, 

T + KF 8(cos T + 2cos 5- 

For the mass of the heaviest particle we then have 



= i u = i ( — + 7 -^-8(cos — + 2 cos -) ) + 0(A 2 ) . (4.5) 



,, o u 
m,2 = lm\ cos — ~ 



2m\ cos 



m 



[!(■ 


Ac 


(m?) 2 


Ac 


7T / 


8 sin — 

5 V 


{m\) 2 



8 sin |fl + 2 cos 2 y ) 



9 2vr\ 
1 + 2 cos 2 1 



(4.6) 



5 / 

where denote the unperturbed value of the masses. By using eq. (C.3) it is easy to 
check that the above expression indeed coincides with the value directly obtained from 
eq. ( p.!6| ), i.e. 

ml = (m°) 2 + 2AF 22 (ivr) . (4.7) 

In addition to this consistency check, a direct test of the theoretical predictions for 
the variations in the spectrum of the theory under the perturbation can be obtained by 
the so-called "truncation method" P5| . The basic idea of this approach is to study the 
theory on an infinitely long strip of width R (the linear spatial volume) with periodic 
boundary conditions. After choosing a Hilbert space basis of eigenvectors of the ultravio- 
let conformal Hamiltonian, all the matrix elements of the perturbed Hamiltonian on this 
basis can be exactly computed. The off-critical spectrum can then be found by numerical 
diagonalisation on a truncated Hilbert space containing a suitable number of states^. 
Since the method does not rely in any way on integrability, it can be applied to any 
perturbation of the conformal point. In particular, for the theory defined by the action 
the energy levels must have the scaling form 

E i (R,g,X) = ^f i (m 1 R; X ) , i = 0,1,2... . (4.8) 



11 All these steps can be performed by means of the algorithm developped in |3£ 
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At very short distance scales, the critical fluctuations are expected to dominate so that 
the spectrum must coincide with that of the conformal point given by |29| 



Ei~^(xi--^ , mi R«l, (4.9) 

where c denotes the central charge and Xi the scaling dimensions of the scaling fields in 
the underlying conformal theory. In the infrared limit, on the other hand, one should 
recover the spectrum of the massive theory on the plane and therefore the energy levels 
are given by 

Ei~ £ vac (g,x)R + mi(g,x) , m l R»l (4.10) 

where the first term takes into account the vacuum bulk energy contribution and m, 
denotes the mass-gap of the i-th leveQ The first energy levels in the spectrum of the 
pure $13 perturbation of the conformal model M( 2 j) obtained by including all conformal 
states up to level five in the Verma modules are shown in Fig. 3. Starting from the 
bottom, the first four parallel lines are easily identified as the vacuum energy level, 
the one-particle energy level of the lowest particle A±, the one-particle energy level of 
the heaviest particle A 2 and the two-particle threshold energy line, respectively. The 
remaining levels will be part of the continuum in the infinite volume limit. The measured 
values of the mass ratio m\/m\ ~ 1.61 and the vacuum energy density £® ac ~ — 0.13(m5) 2 
are in good agreement with the theoretical expectations of the corresponding integrable 
theory. As for other spectra of integrable models [B5|, 137], |38], Rj| , several crossings of the 



energy levels are expected and they are indeed observed for this integrable deformation 

of the M(2,7) model. When a small $1,2 perturbation is switched on (keeping g fixed), the 

first qualitative effect consists to resolve the degeneracy present at the crossing points, 

i.e. all energy levels have now the tendency to repel each other. This feature can be 
12 The above asymptotic expression only applies in the ideal situation where we discard the effects 
induced by the truncation of the Hilbert space. Although one can also take into account these extra 
numerical effects and refine consequently the above formulas, the long practice with integrable massive 
perturbations has shown that as far as the low lying energy levels are concerned, they are not much 
affected by truncation effects so that they can be obtained with a remarkable accuracy on a large range 
of R by just including few conformal states in the calculation. This observation obviously helps to speed 
up the numerical work. 
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interpreted as the breaking of the initial integrability of the model since it is well known 
that, in absence of any higher symmetry, two hamiltonian lines cannot generally cross. 
At a quantitative level, for small values of A, both the separation and the slope of the 
lines vary linearly with it. The first order variations of these quantities are measured to 
be 

- — ~ 3.74 , 

07711 

3£vac r\ r*7 

— ~ 0.67 mV . 

They must be compared with the theoretical predictions for the same quantities obtained 
by using eqs. ( |3.16| ) and ( |3.17| ) 

2 4 cos — sm — 1 = 3.8541 



5m 1 m% fH' 2) (in) V 5 5 



SEvac o <O|01,2)O> COSf ncoo1n o 

— = mV , ' = 2_ = 0.68819.. mV . 

6m, FH' 2 \in) 2sinf 1 

The agreement between the theoretical estimates and the measured values is therefore 
quite satisfactory, the discrepancy being of the same order of the numerical error intro- 
duced by truncation effects. 

The truncation method obviously allows us to study the spectrum of the theory for any 
values of the coupling constants g and A and to explicitly test several theoretical assump- 
tions. We have for instance verified the validity of the scaling law fl3.11|) for the spectrum, 
which in this theory always consists of two distinct one-particle states whose mass ratio 
smoothly interpolates between the values m,2jvn\ = 2cos| and m%jm\ = 2cos^ (rel- 
ative to the limits % — > and \ ~^ °°i respectively), as shown in Figure 4. Obviously 
for large values of \ t ne physical properties of the model can be no longer theoretically 
predicted by the first terms of perturbative series based on the $1,3 deformation. In 
particular, to study the model in the limit x — > 00, it is evident that instead of including 
higher order terms of the perturbative series based on the $13 deformation, it would 
be more convenient to make use of the other perturbation theory based on the second 
integrable operator $ 12 deformation. 
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5 The Scaling Region of the Two— Dimensional Ising 
Model 



Aim of this section is to extract through the Form Factor perturbative techniques some 
information about the scaling region nearby the critical point of the two-dimensional 
Ising model described by 

A = A( 3A )+T J d 2 xe(x) + h J d 2 xa(x) , (5.1) 

and also to study the model by means of the truncation method. As expressed by fl5.1| ), 
the most general off-critical realization of the Ising model is given by a perturbation of the 
simplest unitary minimal model Mi^m with a linear combination of its two only relevant 
primary fields, i.e. the energy density $ 13 (x) = e(x) of scaling dimension x £ = 1 and the 
magnetization &i^(x) = cr(x) of scaling dimension x a = 1/8. The conjugated couplings 
t and h are respectively interpreted as the deviation from the critical temperature T c 
and a constant magnetic field, with physical dimensions given by r ~ M and h ~ M 15 / 8 , 
where M is a mass scale. The Zamolodchikov's counting argument f2| ensures that the 
theory (|5.1|) becomes integrable as far as one of the two coupling constants r or h is set 
equal to zero. The integrability of the purely thermal perturbation has been known for 
long time and can be neatly reformulated in terms of field theory of Majorana fermions 



[[12], [15], |22], [23], [24]]. The integrability of the magnetic deformation has been established 
by Zamolodchikov || in the recent past and, for the original features of his findings this 
result may be regarded as a direct success of the application of field theoretical methods 
to statistical mechanics. 

The action ( p.l| ) defines a one parameter family of field theories which can be labelled 
by the dimensionless combination x = t|/i| _8//15 G (— oo,+oo). Previous investigations 
have shown that the particle content of the model must undergo drastic chang es as a 



function of x 2"6fl - As we will see in the following, this prediction can be directly 
confirmed by the numerical determination of the spectrum for different values of x ( see 
Figures 7). The reason for such behaviour of the spectrum will be easily understood after 
reviewing some exact results obtained for the two integrable directions. 



25 



5.1 Thermal Deformation 

Consider first the case h = 0. In the high temperature phase, i.e. % = +°°) the theory 
has an unique vacuum state and the spectrum consists of a single particle A of mass 
m = 2ttt and multi-particle states theoreof. This is clearly confirmed by the truncation 
method calculation shown in Figure (7. a): starting from the bottom, the sequence of 
energy levels is given by the ground state, one-particle level and then all the multi- 
particle lines of the continuum, the lowest of which corresponding to the two-particle 
threshold. The massive excitation present in this phase of the model can be equivalently 
regarded as a free neutral fermion or as an interacting boson (with scattering amplitude 
S = — 1) created by the magnetization operator. Invariance under the spin reversal 
implies that the "order" ("disorder") field a(x) (/i(x)) couples to the states with an odd 
(even) number of particles only. The Form Factors for the two operators are given by 
the unique expression fl2|, [15|1 

n n n 

F {n) {e u ...,e n ) = {Q\0{Q)\a{e 1 )...a{e n )) = (#~ X J ntanh^— -i , (5.2) 

i<j 1 

with n odd for 0(x) = cr(x) and even for 0(x) = /z(x) and [a] denoting the integer part 
of a. 

The situation in the low temperature phase, i.e. x — —oo, is easily deduced from 
the previous one by duality. In fact, the excitations are now created by the disorder 
operator fi(x) and must be interpreted as kinks and anti-kinks of mass m interpolating 
between the two degenerate vacua coming from the spontaneous symmetry breaking of 
the Z2 invariance of the model. The Form Factors of <j{x) coincide with those of the 
operator /i(x) computed in the high-temperature phase T > T c . Hence, the magneti- 
zation operator o(x) now couples only to states with even number of excitations. The 
numerical determination of the spectrum is shown in Figure (7.g). The lowest two lines 
are clearly the degenerate ground state levels. They approach each other exponentially 
in the crossover region, AE ~ exp(— mR), where m is the mass of the kink responsi- 
ble of the finite-size tunnelling effect between the two vacua. Since periodic boundary 
conditions were adopted on the strip of width R in the numerical determination of the 
spectrum, in Figure (7.g) there are no energy levels corresponding to odd number of kink 
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states. In particular, the energy line associated to the single kink state is absent and 
therefore the lowest energy level above the ground state energies is given in this case by 
the kink-antikink threshold line. 

5.2 Magnetic Deformation 

Far richer is the situation for the other integrable direction obtained when r = and 
h 7^ 0, i.e. x = 0- m fact, Zamolodchikov has shown that in this case the spectrum 
consists of eight different species of self-conjugated particles A iy i = 1, . . . , 8 with masses^ 
i 



mi 


= M(h), 






m 2 


71 

= 2mi cos — = ( 
5 


1.6180..) 


mi , 


m 3 


7T 

= 2777! COS — = 


(1.9890.. 


) mi , 


m,4 


= 2777 2 COS = 

30 


(2.4048.. 


.) 777-1 , 


m 5 


2tt 

= 27770 COS = 

15 


(2.9562. 


.) 777,1 , 


m 6 


7T 

= 27770 COS = 

30 


(3.2183.. 


) mi , 




= 47772 COS — COS 
5 


^ = (3. 
30 v 


8911..) 777i 


777 8 


= 4777o cos — COS 

5 


— = (4.7833..) mi 
15 V ; 



(5.3) 



Notice that only the first three particles of the above spectrum lie below the lowest thresh- 
old at 2t?7i. The remaining particles are however stable since integrability prevents the 
possibility of inelastic processes, in particular decay processes. The S-matrix describing 

the interactions between the eight particles was exactly determined in 0. These results 
13 Within the standard CFT normalization of the magnetization operator, obtained by the equation 
< a(x)a(Q) >= — ^-r (\x\ — > 0), the overall mass scale M(h) has been exactly determined in Jlo| , 

\x I ^ 

M(h) = Chrs where 

4sin f r(i) / 4 ^r(f)r»(g) 
r(|)r(A) ^ r(i)P(^) 



4.40490858. 
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were used in refs. [1 1] to implement the bootstrap equations for the Form Factors of the 



model and to compute in particular the one and two-particle matrix elements of the op- 
erator a(x). As expected from the fact that neither the action nor the S-matrix of the 
model exhibit any internal symmetry, the magnetization operator was found to couple to 
all the eight particles in the spectrum, namely F? = (0\cr\Ai) ^ for i — 1, ... ,8. 

The finite size spectrum of the Ising model in a magnetic field at T = T c as obtained 
with the truncation method is shown in Figure (7.d). Starting from the bottom, five 
parallel lines are clearly identified in the large volume limit. They correspond to the 
ground state, the three lightest particles and the first two-particle threshold, respectively. 
The remaining single particle states are of more difficult identification since they are 
placed among the higher threshold lines. Notice that they cross in several points the 
"momentum lines" (i.e. lines corresponding to states containing two or more particles 
not at rest with respect to each other) converging toward the threshold in the infinite 
volume limit. Level crossing in the energy spectrum is always related to the presence of 
symmetries and in the present case strongly supports the integrability of the theory. 

5.3 McCoy— Wu Scenario 

After having discussed individually each deformation, let us consider the more general 
case described by the action (|5.1|). Following McCoy and Wu P5| , let us introduce the 
Fourier transform of the spin-spin correlation function 

G(p 2 , X )= Jd 2 xe^{a(x)a(0)) . (5.4) 

According to the discussion of the previous subsections, the leading infrared singularities 
of this function in the complex p-plane for \ = +oo, 0, — oo are those shown in Figs. (6. a), 
(6.d) and (6.g), respectively. Since these three situations correspond to different values 
of the couplings in the same QFT defined by the action (|5.1|), it must be possible to 



interpolate continously between them by moving along the semicircular path drawn in 
Figure 5. To describe the change of physical properties of the theory along the path, 
the following scenario has been proposed |25|] : starting from the point a in Figure 5, 



the spin reversal symmetry will be broken as far as a magnetic field h is turned on 
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and consequently, even particle thresholds immediately appear, as shown in Fig. 6.b. 
Moreover moving toward the point d, additional poles emerge from unphysical sheets of 
the Riemann surface through the branch cuts and become physical bound state poles. 
From Zamolodchikov's solution, we know that at the point d there are precisely three 
such poles below the first two-particle threshold. Their number however continues to 
increase when moving toward the negative r axis until altogether they coalesce at point 
g and give rise to the branch cut starting at 2im in Fig. (6.g). The lowest lines of the 
numerical spectra corresponding to the values [] X = +00, 1.32, 0.16, 0, —0.16, —1.32, 
— 00 are shown in the sequence of Figures 7 and we will comment on them below. 

The coalescence of the poles is quite a striking phenomenon. It was quantitatively 



discussed in ref. |25| where the knowledge of the n-point correlators for the purely thermal 
Ising model was exploited to study the behaviour at small h of the spin-spin correlation 
function expressed in the form 

00 1 f 

(cr(x)a(0)) conrii h = J2 ~ hn / • • • d 2 x n (cr(x)a(0)a(xi) . . . cr(x„)) co?m> ^=0 ■ (5.5) 

It was found that when the magnetic field is switched on, the two-particle branch cut of 
Fig. (6.g) breaks up in a sequence of poles located at 

(2 + h 2/3 -f 2 k /3 )im , (5.6) 

where 7^ are the positive solutions of 

J i (b k ) + J -i (b k ) = • 

and J v {x) is the Bessel function of order v. The above values coincide with the eigenvalues 
of the Schroedinger equation for a particle in a central linear potential. The physical origin 
of this result is that the presence of a small magnetic field breaks the degeneracy of the 
two vacua in the low-temperature phase inducing a linear confining potential between 
the kinks0. It is clearly evident then that the magnetic field induces a drastic change in 

the structure of the low-temperature phase of the model: in fact, as soon as the magnetic 
14 The quoted values of x refer to the standard CFT normalisation of the fields e(x) and a(x). For a 

generic field tp of scaling dimension x, this is defined by the condition (cp(r)(p(0)) — > r~ 2x , r — ► 0. 
15 One can easily convince himself of this feature by using a semi-classical analysis of the kink states. 
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field is switched on, the kink configurations are no longer asymptotic states of the field 
theory and the original kink-antikink pairs collapse into a sequence of bound states. This 
happens for any value of the magnetic field h, however small. The lightest of these bound 
states has a mass equal to 2m + 0(h 2 ^ 3 ), so that the first branch point is now located at 
4im + 0{h 2 / 3 ). 

5.4 Form Factor Perturbation Theory 

Let us now turn to the discussion of the Ising model within the perturbative frame- 
work considered in this paper Q. Our main interest will be to develop the perturbative 
picture around the integrable magnetic direction (x = 0), for which no other approach 
is presently avalaible. Before doing that, we will however briefly comment on the case 
h ~ 0, distinguishing the two situations T > T c and T < T c . 

At T > T c , the magnetization operator er(x) couples to the odd-particle states only. 
Hence, the opening of inelastic channels (e.g. the production process AA — > AAA) with 
the consequent breakdown of integrability can be easily checked at first order in h but 
the correction to the mass or the variation of the elastic scattering amplitudes occurs 
instead only at the second order in h. The breaking of the integrability implies that 
the numerical spectrum no longer presents in this case crossing of the energy levels, a 
feature which, although not always clearly visible in Figs. 7, is however confirmed by the 
numerical data. 

In the low-temperature phase T < T c the spectrum of the theory is known to undergo 
a qualitative and non-analytic change when the magnetic field is switched on. Although 
these considerations naturally suggest the failure of any finite order of a perturbative 
approach to study this phenomenon, it is nevertheless instructive to see how such failure 
manifests itself in our method. To this aim, it is sufficient to attempt the computation of 
the correction to the mass of the kink, an excitation which we know no longer exists as an 
asymptotic state in the perturbed theory. According to the general formula ( |3.16| ), the 



first order correction is proportional to the two-particle form factor of the magnetization 



16 



For other perturbative studies of the Ising model field theory see refs. M, Ell 
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operator computed for a rapidity difference equal to in. Since this form factor is equal to 
tanh ^ 2 (see eq. (|5.2| )), it has a pole at 9i — 9 2 = in and therefore an unbounded correc- 
tion is obtained for the mass of the kink. The general theory of Form Factor predicts that 
such "kinematical" pole is present in the two-particle Form Factor (O\O(O)\A(6i)A(02)) 
only if the operator 0(x) is non-local with respect to the field which creates the particle 
A (as it is the case for the fields cr(x) and fi(x) in the Ising model). This circumstance 
suggests therefore that a confinement phenomenon of the type described above for the 
Ising model has to be expected each time that the perturbing operator and the fields 
which create the particle excitations in the unperturbed theory are not mutually local. 
It would be interesting to check this prediction in other statistical models where the 
integrable spectrum is given by kink excitations. 

In the case of the Ising model, the numerical determination of the spectrum confirms 
the above scenario. From Figure (7.f), we see in fact that the first effect of h consists 
in the removal of the degeneracy of the lowest as well as of the higher eigenvalues. The 
two originally degenerate ground state levels have now been splitted into: (i) a unique 
ground state energy line and (ii) an excited state whose mass gap diverges in the large 
volume limit R — > oo. The degeneracy of the lowest threshold line has also been lifted 
giving rise to a sequence of one-particle energy levels, the lowest values of which can be 
checked to be in reasonable agreement with eq. ( |5.6| ). By increasing the magnetic field h, 
the above characteristic feature of the spectrum are futher enhanced, as seen by following 
Figs. (7.e) and (7.f) in reverse order: the divergent energy line of the initial degenerate 
ground states meets all other lines at smaller values of R (and therefore it decouples 
faster from the remaining spectrum) whereas the other eigenvalues start to assume the 
structure and the values predicted by the Zamolodchikov solution for T = T c , reproduced 
in Figure (7.d). Observe that, while the eigenvalues have the typical repulsive behaviour 
of a non-integrable situation all along the path from the point / to e, on the contrary they 
cross each other once we have reached the integrable situation of the pure magnetic axis. 
Moreover, the fact that in the pure magnetic case the approach to the asymptotic value 
of the masses is reached from below can be simply interpreted as the reminiscence of the 
divergent line coming from the degenerate ground state energies of the low-temperature 



31 



phase. 

Let us now analyse more closely and in quantitative terms the field theory defined in 
the vicinity of the magnetic axis. Perturbation theory around the magnetic axis (T = T c 
with h 7^ 0) obviously requires the knowledge of the Form Factors of the energy operator 
e(x) in the Zamolodchikov field theory involving the eight massive states (i — 1, • • • , 8) 
associated to the pure magnetic model. These Form Factors can be exactly computed by 



using the same bootstrap approach used in ref. |TT] for the magnetization operator. The 



detailed discussion on the computations of these quantities will be presented elsewhere 



[pl| and here we only quote the results needed for our present purposes 



(0|e(0)|0) = m 1 , 

F^iir) = (O|e(O)|Ai(0 + i7r)Ai(0)) = -17.8933. .mi , 
F 2 £ 2 (ivr) = (0\e(0)\ A 2 (6 + i7r)A 2 (6)} = -24.9467..mi , 
FJgOr) = (0\E(0)\A 3 (9 + irr)A 3 (9)) = -53. 6799. .mi , 
F 4 £ 4 (» = (O|£-(O)|A 4 (0 + ztt)A 4 (0)) = -49.3169. .mi . (5.7) 

The first line of the above relationships should be meant as expressing the normalisation 
of the energy operator e(x). When the energy perturbation is switched on, the spectrum 
in Fig. (7.d) undergoes continous deformations which can be followed by the truncation 
method. In particular, for |r| ~ \T — T c \ sufficiently small and h fixed, the variations 
in the energy levels are linear in r and their direct measurement with the truncation 
method is given by 



Smi 
Sm,2 
5m i 
5m 3 



-0.05 m? 



0.87 , (5.8) 
1.50 . 



5 m i 

The theoretical predictions for the same quantities are obtained plugging the values 
into the first order formulae (|3.17 ) and (|3.16j ) 



$£yqc (OHO) o nn^o o 

— = - ., r my = — 0.0558. .mV 

5mi Ff^iir) 
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— ^ = 22 L —L = 0.8616.. , 5.9 
dmi Fii^m) 

pL = ?kM 4 = 1.5082.. . 

The agreement between the theoretical and numerical estimates should be regarded as 
quite satisfactory. 

More dramatic is the effect of the perturbation on the five particles above the thresh- 
old. We have already mentioned that their stability at r = must be considered a 
consequence of the integrability of the theory which prevents any kind of inelastic pro- 
cess. They are then expected to decay in the perturbed, non-integrable theory. This is 
easily seen to be the case by analysing the effect of the perturbation on the finite size 
spectrum of the model. Indeed, since integrability is lost under the perturbation, level 
crossing in the spectrum is no longer allowed at r / 0. This means that the energy 
levels are now forced to repel each other at the former crossing points. Since each line 
associated to a particle above threshold crossed an infinite number of momentum lines 
at r = (Fig. (8. a)), it immediately becomes a "broken line" (almost) parallel to the 
threshold when the perturbation is switched on (Fig. (8.b)). This is exactly the signature 
of unstable particles in the finite volume; more precisely, the difference in the slopes of 
the broken line and the threshold is proportional to the width of the resonance |34| . 



In perturbation theory, the decay of the particles above threshold manifests itself 
through the appearence of a negative imaginary part in the mass. This is an effect which 
occurs though at the second order in r. To be specific, consider the lightest unstable 
particle A±. It follows from (|3.16| ) and (|5.7f) that the first order correction is real and is 



given by 

^ = ^4 = 1.1460... (5.10) 

The computation of the second order correction requires the sum over a complete set of 
intermediate n-particle states and involves an energy denominator. Expressing the latter 
in the usual form 

g -i„ +i . = p (^;)-' rf(£ - £ " ) - (5 - n) 

one immediately concludes that the principal part gives rise to the second order correction 
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to the real part of the mass while the delta function originates an imaginary part receiving 
a contribution only from the intermediate states whose total energy-momentum equals 
that of the external particle. For the particle A A we have 

-2 

2m1m1 sinh 9* 



Imml ~ - ~ Q l.^M M*)m\MnM-0*M , (5-12) 



where 6* = \circcosh^\ 



6 Conclusions 

In this paper, we have shown that once the complete dynamics of two-dimensional exactly 
solvable models is known, this can be also extremely useful for investigating the structure 
of the quantum field theories close to the integrable directions. In particular, we have 
derived a perturbation theory for the non-integrable QFT based on the Form Factors of 
the exactly solvable relativistic models and we have considered the lowest term of the 
series in order to study the effects of a small perturbation which breaks the integrability 
of the unperturbed theory. We have discussed the variation of the mass spectrum, the 
shift in the vacuum energy density and the correction to the elastic part of the S'-matrix, 
and we have successfully checked them against their numerical estimates obtained by the 
truncation method in two models with an underlying CFT, namely the minimal model 
M( 2> 7) and the Ising model. 

It is worth pointing out that, specialising the perturbative formulas to the trivial case 
of a parallel perturbation along the original integrable model, one may obtain as a by- 
product an infinite set of useful identities for the Form Factors of the stress-energy tensor 
of the original integrable field theory. While this appears to be as a convenient method 
of deriving them, at the same time it explains the reason of their validity because they 
can be regarded as consistency equations for the perturbative scheme built up on the 
Form Factors. Known features of the Form Factors of the integrable relativistic models 
are then deeply inter-related with the general structure of quantum field theories. 

The methods illustrated in this paper may apply of course to other interesting physical 
situations in addition to those considered here. We would like to mention, for instance, 
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two non-integrable lagrangian field theories which would be interesting to analyse in 
terms of their Form Factor perturbative series. The first one is given by 



1 /a ^2 m2 



1 - - J {2e 9v + e" 2 ^) + \{e 9V + e _w ) 



By varying the parameter A in the interval [0, 3], we can interpolate between the Bullogh- 
Dodd and the Sinh-Gordon models. Both these field theories are separately integrable and 
the matrix elements of their local operators have been computed in |T7| , 18 1. Notice that 
the Z 2 symmetry of the Sinh-Gordon model is always broken along the entire interpolating 
trajectory and it is only recovered at the end point A = 3. 

The second non-integrable lagrangian model is obtained by adding higher "harmon- 
ics" to the original Sinh-Gordon interaction, the simplest example of this class of non- 
integrable models being the so-called double Sinh-Gordon model 

L = -(d^) 2 - — [(l-4A)cosh^ + Acosh2^] . (6.2) 

By varying A in the interval [0, ~], we can interpolate in this case between the integrable 
field theories given by the Sinh-Gordon theory with coupling constant g and 2g. 

Notice that for both examples ( |6.1| ) and ( |6.2| ), the integrability of the two theories at 
the two extremes of the interpolation interval implies strong constraints for the ordinary 
perturbation theory based on the coupling constant g. It would be then interesting to 
learn more about those models by comparing the results obtained in ordinary Feynman 
diagram perturbation theory with those derived by the perturbative series based on the 
Form Factors. 
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Appendix A 



A simple example of quantum mechanics will help to clarify the simultaneous relevance 
of all the inelastic thresholds in any energy interval. Consider a one-dimensional inelastic 
scattering process of a particle of mass m which hits a target with n + 1 internal states 
of increasing energies (i = 0,1, ... ,n), where the excited states | i > (i = 1, . . . , n) 
are given by | i >= a\ | > (aj and denote the creation and annihilation operators, 
with standard commutation relations). The states of this system (target plus particle) 
are described by the wave function | \I> >= YJi=o a \ I > 4>i(x) and as Hamiltonian of 



the system we choose 

n r 2 n n 

(A.l) 



p 2 n h 2 



1m ^ 1m 

Z " L i=0 Z " L L i=0 i=l 

U rules the elastic transition amplitudes | i >^\ i > whereas Ui are related to the 
inelastic reactions | ><->| i >■ For this toy- model, it is quite simple to determine the 
phase-shift 5 and therefore the S'-matrix element S = e 2tSo relative to the elastic channel 
| >— >| >. This is given by the formula 

tan<5 = ^-V--^ — , (A.2) 

where 

The phase shift is rea l when the energy E is below the first threshold E\ and complex 
above but it is worth to note that also in the elastic region < E < E\, its value is 
determined by all the inelastic parameters of the problems, i.e. U\ and k\. 



Appendix B 

As it is well known, a very effective algebraic description of integrable theories can be ob- 
tained in terms of the Faddeev-Zamolodchikov (FZ) creation and annihilation operators, 
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Z+(9) and Z a {6). They satisfy the algebra 

z:(e 1 )z+(e2) = sSft - e 2 ) z+(9 2 )z^(9 1 ) 

z a {o x )z b {o 2 ) = sSfr - o 2 ) z d {e 2 )z c (e{) (b.i) 

z a (Q x )zt{Q 2 ) = s£{e 2 - e 1 ) zi{d 2 )z c (Q x ) + 2ns ab s(e 1 - e 2 ) 

which can be regarded as a generalisation of the canonical commutation relations. The 
operator Z a {6) annihilates the vacuum while Z+{9) creates a particle of tipe a with 
rapidity 9. The space of states is generated by 

|oi(eO---an(e„)> = z+fa) . . . z+(e n )\o) , 

{a n {e n )...a 1 {e 1 )\ = (o\z an (e n )...z ai (e 1 ) . (b.2) 

The physical asymptotic states can be selected through the following ordering prescription 
over rapidities: the states ( |B.2| ) are "In" states if 9\ > 9 2 > . . . > 9 n , and "Out" states if 

e 1 < e 2 < . . . < o n . 

It is easily checked that the unitarity and Yang-Baxter equations for the S-matrix 
can be obtained requiring respectively the consistency under double application and the 
associativity of the FZ algebra (|B.1|) . Concerning the remaining fundamental property 
of the S-matrix, namely the crossing relation 

Sf b (0) = St(i7r - 6) , (B.3) 
let's formally perform in the last equation in (|B.1|) the analytic continuation 9 2 — > 9 2 +m. 



Then, substituting particles b and d with their antiparticles b and d and using eq. (B.3|) , 
we obtain 

Z a (e 1 )Z- b + {9 2 + m) = Sf d (9) Z±(9 2 + m)Z c {9 x ) , [9 X + 9 2 ) . (B.4) 



Comparison with the second equation in ( B.I ) suggests the identification 



Z l + {9 + m) = -Z b {9) . (B.5) 

The proportionality constant —1 has been chosen in order to fit the canonical cases 
corresponding to S^{9) = ±1. Indeed we recall that the Fourier decomposition of a free 
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(e.g. bosonic) field 



1 r AU 1 r 

V (x) = -j — [a(k)e-^ + a + (ky k * 



immediately leads to the representation 



a(k) = -- j dx 1 e ikx d ip(x) = f(k) , 
a+(k) = lJdx 1 e- ikx d <p(x) = -f(-k) 



Consider the two matrix elements 



(B.6) 



coFZ^K • • • A) = (o\[Zco(e), 0(0)]!^), • • • , a n (e n ))+(o\o(o)z co (e)\ai(0il •••>««( 

(B.7) 

F£ ai ... an (0, Oi,...,0 n ) = (0|[O(0), Z+W^ijh), . . . , a n {9 n )) . (B.8) 

We can use the identification ( B.5 ) and the FZ algebra to recognize that the first matrix 
element can be written as 



F au..a n (^U • • • A 



CO 



8=1 



F? (t 

1 COOL. .a„ V 1 

k=l 



m,9 x , ...,9 n ) + 



(B.9) 



F dl ... di _ iai+1 ...a n {9l, ■ ■ ■ ^i-U • • • > On) 



This result is represented pictorially in fig.Al. As an illustration of this crossing procedure 
we write down esplicitely the two matrix elements involved in the evaluation of the mass 
and amplitude variations ( |3.16| ) and ( |3.21| ) 



h F°{6 2 \e{) = F&(6 2 + m, 9 X ) + 2n5 ab 5(9 1 - 2 )(O\O\O) 



(B.10) 



aFg{9 z , 9±\9 U 9 2 ) = Fg ab {9z + in, 9, + in, 9 U 9 2 ) + 
27r5(0i - 9 3 )Sf a (in + 9 4 - 9 3 )F e ° b (9 4 + in, 9 2 ) 
+2ir8{9 2 - 9 A )S e a d b (9 1 - 9 2 )Fg{9 3 + in, 9,) + 



+(2n) 2 5(9 2 - 9,)5(0i ~ 3 )S, 



9 2 ){0\O\0) 



(B.11) 



In the last equation we excluded the "resonances" # 4 = 9\ and 9 3 = 9 2 . 
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Appendix C 

In this section, we collect the needed FF of the operator $12 relative to the $1,3 defor- 
mation of the minimal non-unitarity model M^,7)- They are given by 

< 0| *- |0 > = 2^f ' (ai) 

5 

^) = -c 4sin-^, (C.2) 

2 7 r (coshg + 4cosfsin 2 f)F— (0) 

F 22 {9) - -c 8sm- Vl{e)Vl{e)V±{e) > ( C - 3 ) 



and 



oi n ^ — ^ Tpmin (fj \ 



In the above formulas, the functions 



and 



where 



F™ n (9) = -ismh^Giie) (C.5) 
2 5 



F™ n (0) = -ismh-G2 T (9)G±(9)G±(9) , (C.6) 

^5 5 5 



G a (#) = exp 



+°° cosh(a — |)i . 2 (in - 9)t 
t cosh I sinh t 2tt 



(C.7) 



solve the monodromy problem for the form factors, while the functions 

appearing in the denominators are in coincidence with the bound state structure of the 
theory. The function Q in the four particle form factor is defined as 



cosh 9ik + 4 sin' 



2tt 



, #12 , #34 
cosh — — cosh — — 
2 2 



7T 

cos — 
5 

COsh ^(0 31 + 32 + 34 ) + COsh i (041 + #42 + M 



, #13 , #24 

cosh — — cosh — — 
2 2 



014 , #23 

cosh — — cosh — — 
2 2 



+ 



COsh i(0 12 + 013 + 014) + COsh -(0 21 + 23 + 24 ) + 
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Notice that the constant Cq in front of all Form Factors depends on the (non-universal) 
normalisationF] of the field $12 (x). We can get around the difficult problem of determin- 
ing this constant by considering universal ratios. 



17 The value of cq corresponding to the standard CFT normalisation, identified by the two-point 
function behaviour < < &i2(^c) < &i2(0) > — -> x~ 4Al2 , can be hardly determined without resumming the 
entire spectral series of the two-point function. This turns out to be in general quite a difficult problem. 
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Figure Caption 

Figure 1 . Analytic structure of two-body S'-matrix in non-integrable field theories. 
The crosses indicate the location of the bound state poles whereas the thick lines 
the branch cuts originating at the thresholds. 

Figure 2 . Renormalization Group trajectories associated to the field theories: inte- 
grable directions (continous lines) and non-integrable one (dashed line). 

Figure 3 . First energy levels of the $1,3 integrable deformation of the conformal model 
M( 2 ,7) as functions of the strip width R. 

Figure 4 . Mass ratio of the non-integrable deformations of the conformal model M^,7) 
as function of x- 

Figure 5 . Phase space in the vicinity of the critical point of the 2-d Ising model. 

Figure 6 . Analytic structure in momentum space of the spin-spin correlation function 
of the Ising model relative to the points (a, . . . , g) in Fig. 5. 

Figure 7 . Numerical spectrum of the Ising model relative to the points (a,...,g) in 
Figure 5. 

Figure 8 . (a) The energy level corresponding to a stable particle above threshold crosses 
several momentum lines in the integrable theory; (b) removal of the crossing points 
implying the decay of the particle above theshold in the perturbed, non-integrable 
theory. 

Figure 9 . Crossing relation for the Form Factors. 
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